QED effects in molecules: test on rotational quantum states of H 2 
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Quantum electro-dynamic effects have been systematically tested in the progression of rotational 
quantum states in the X Jit ,v = vibronic ground state of molecular hydrogen. High-precision 
Doppler-free spectroscopy of the EF 1 ^ — X 1 E^" (0,0) band was performed with 0.005 cm -1 
accuracy on rotationally-hot Hb (with rotational quantum states J up to 16). QED and relativistic 
contributions to rotational level energies as high as 0.13 cm -1 are extracted, and are in perfect 
agreement with recent calculations of QED and high-order relativistic effects for the H2 ground 
state. 

PACS numbers: 31.30.J-, 33.20.Lg, 33.20.Sn 



Quantum electrodynamics (QED) has been hailed as 
the most successful theory in physics as its predictions 
are in remarkable agreement with a variety of extremely 
precise experiments. Precision tests of QED include free- 
particle systems (e.g., the anomalous magnetic moment 
of the electron [l| and of the muon Q) and bound atomic 
systems, such as the (one-electron) H atom [3[, heavy hy- 
drogenlike ions (e.g., U 91+ [4J), and the (two-electron) He 
atom [5j. Recent reviews of QED theory and precision 
tests in simple atoms and ions can be found in Refs. [6J,lZ|. 
Progress has been made in QED calculations for the 
smallest (one-electron) molecular ion H 2 + Q, and sys- 
tematic high-resolution spectroscopic investigations have 
been proposed based on cooled ions in a trap 9]. Re- 
cently, the problem of QED calculations (also including 
high-order relativistic corrections) in the smallest neutral 
molecule H2, the benchmark system of molecular physics, 
has been addressed by Pachucki and co-workers. The ex- 
perimental determinations for the ionization and dissoci- 
ation energies of molecular hydrogen and its isotopomers 
10~ 8 level of accuracy for H 2 [lfj], D 2 ll| and 



to the 
HD 



12L Il3f had stimulated calculations of these quan- 



tities. The ab initio calculations of the FEand D 2 dis- 
sociation energies 



141, and that for HD 



15| . for which 

relativistic and QED effects for the lowest J = level 
were computed, are in excellent agreement with exper- 
imental results for the dissociation energies. The QED 
calculations were thereupon extended using the theoret- 
ical framework in Refs. [14J, [l5| to the full set of rovi- 
brational levels in the ground states of H 2 , D 2 and HD. 
When combined with the updated nonrelativistic ab ini- 
tio calculations including adiabatic and nonadiabatic ef- 
fects [16[ , a full set of QED rovibrational energies for H 2 , 
D 2 , and HD are now available 
In the present investigation 



laua 



we pursue a systematic 



study of QED effects (including higher-order relativistic 
and radiative effects) in a progression of 16 rotational 
quantum states in the X 1 ^, v — ground state of H 2 . 
For this purpose, a nonthermal H 2 population distribu- 
tion was produced via ultraviolet(UV)-induced photodis- 
sociation of hydrogen bromide (HBr) and a subsequent 
chemical reaction: H + HBr -> Br + H 2 (J) [lq . thus 
forming rotationally hot H 2 (J) with quantum states de- 
tected up to J = 16. Here we follow similar procedures 
employed by Heck et al. [19( for the production and (low- 
resolution) spectroscopic studies of hot D 2 . 

The principle behind the derivation of the accurate 
rotational level energies is based on the precise laser 
spectroscopic measurements of EF 1 ^ — X x 5]+ ,Q(J) 
two-photon transitions (Q lines denote AJ = tran- 
sitions [20|) in the range J = 6 — 16, combined with 
the existing information on the excited states. This ap- 
proach depends on the availability of accurate level en- 
ergies of EF l Y,+, (v = 0, J = 2 - 12) lrpm_studies de- 
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tcrmining EF, J = 0, 1 anchor levels [20|, |21| and per- 
forming Fourier-transform spectroscopic studies on the 
manifolds of excited states in H 2 [22j . By measuring ad- 
ditional O-branch (AJ = -2) and S-branch (AJ = +2) 
two-photon transitions, X 1 E+ ground state energy split- 
tings were derived extending the rotational sequence up 
to J = 16. Experimental relativistic and QED correc- 
tions to the ground state rotational levels are finally ob- 
tained by subtracting the calculated nonrelativistic ener- 
gies, taken from the recent ab initio study of Pachucki 
and Komasa [16j |. 

In the experiment, a pulsed-dye-amplifier (PDA) laser 
system, which is injection seeded by a continuous-wave 
(cw) ring dye laser [23]], is operated at wavelengths be- 
tween 610 — 645 nm. Pulsed ultraviolet radiation of nar- 
row bandwidth in the 203 — 215 nm wavelength range 
is obtained via two-stage third-harmonic up-conversion 
of the PDA-output radiation in nonlinear crystals. The 
method of 2+1 resonance-enhanced multiphoton ioniza- 
tion (REMPI) is used to probe the EF — X transi- 
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TABLE I. Frequencies (in cm ) of two-photon transitions 



in the H 2 EF 1 ^ 
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-0.1 0.1 

Energy + 95 875.406 (cm" 1 ) 

FIG. 1. (Color online) Recording of the EF 1 ZJ-X X E+ (0,0) 
two-photon 0(9) transition. Relative frequency markers are 
provided by the fringes of a stabilized etalon (FSR=148.96 
MHz), while absolute calibration is achieved using the known 
position [25|] of an I2 hyperfine component (marked by *) as 
the absolute frequency standard. 



tions [20|. We note that the production of hot H2 and 
the spectroscopy is performed by the same UV laser pulse 
and occurs within 5 ns. Ions produced via REMPI are 
accelerated through a time-of-flight mass selector and de- 
tected by an imaging system composed of a multichannel 
plate, phosphor screen and photomultiplier, and finally 
digitally registered. 

To obtain Doppler-free spectra, two counterpropagat- 
ing UV laser beams are crossed with the molecular beam 
of HBr, in which rotationally hot H2 is produced. The 
counterpropagating UV beams are aligned using the 
Sagnac geometry [24j to avoid residual Doppler shifts, 
which we estimate to be less than 1 MHz for the in- 
vestigated H2 transitions. The application of a separate 
355-nm ionization laser (pulse delayed by 30 ns with re- 
spect to the excitation laser) enables us to reduce the 
UV radiation intensities, in order to produce the narrow- 
est line profiles (350 MHz) and minimize ac Stark shifts. 
Typically, only 50 — 200 /iJ is used in the experiment 
(for hot-H2 production and two-photon excitation) out 
of the ~ 1 mJ UV radiation produced. For most (strong) 
transitions, we recorded spectra at different UV intensi- 
ties for extrapolating to zero-intensity frequencies within 
~ 50 MHz. dc Stark shifts were found to be negligible. 

To provide a relative frequency scale, part of the cw- 
output of the ring dye laser is directed to an actively- 
stabilized reference etalon with free spectral range (FSR) 
of 148.96 MHz. Another part of the ring laser cw-output 
is utilized in a Doppler-free molecular iodine (I2) satu- 
ration spectroscopy setup for absolute frequency calibra- 
tion [251 ] . A typical measurement result is depicted in 
Fig. [TJ which shows the simultaneous recording of the 
H2 EF — X, (0, 0), 0(9) transition, the stabilized etalon 
fringes, and the hyperfine-resolved I2 calibration line 



estimated to be 0.005 cm 



V X E+ (0,0) band. 



The uncertainty is 
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[i?(82) (10 - 4) ai component at 479 045193.7(5) MHz]. 
Note that in Fig. [1] the energy scale (lower axis) of the 
H2 resonance is exactly sixfold that of the etalon markers 
and the I2 reference spectrum (upper axis). The uncer- 
tainty contribution to the EF — X transition frequen- 
cies from the scan (non) linearity and the I2 calibration 
are estimated to be 5 MHz, respectively. For typical H2 
transitions recorded, the signal-to-noise (SNR) ratio is 
sufficient to obtain the line position to within 10 MHz. 
The largest contribution to the experimental uncertainty 
is systematic, and derives from the estimated frequency 
offset of ~ 150 MHz associated with frequency chirp dur- 
ing the pulse amplification process in the PDA [26j. In 
total, we estimate the experimental uncertainty of the 
transition energies to be 160 MHz or 0.005 cm -1 . 

The measured transition frequencies of Q-, O- and S- 
branch lines are listed in Table |U To validate the accu- 
racy of the present measurements, we have remeasured 
the Q(5) transition energy and obtained agreement to 
within 0.000 35 cm" 1 of the previously measured value 
at higher accuracy (0.00015 cm -1 ) 2_l|. This indicates 
that for this particular frequency range, the chirp effects 
in the PDA are much less than estimated. Nevertheless, 
we retain the conservative estimate (0.005 cm" 1 ) in view 
of deviations found in other wavelength ranges [22J . 

Ground state rotational level energies X(J) are ob- 
tained by subtracting the experimental EF—X transition 
energies from the EF, v — 0, J = 2 — 12 level energies 
as determined by Bailly et al. [22J. For J = 13 — 16, the 
ground state rotational level energies are derived from 
the Q- and O- or S'-branch transition frequencies. The 
combination differences between Q(J) and S(J — 2) and 
between Q{ J) and 0(J + 2) yield ground state energy 
splittings that validate the assignments of the transitions 
and the EF(J) level energies of Bailly et al. [22J; more- 
over, these measurements demonstrate that EF(13) was 
missassigned in Ref. [22j. The resulting values are listed 
in Table HI! which also include level energies X(J = 2 — 5) 
derived from Refs. [20j, |21| and the X(J = 1) level en- 



TABLE II. Rotational energies X(J) of H 2 X, v = levels 
with respect to the X(J — 0). The relativistic and QED 
corrections X re \+QST>(J), with respect to X ic i+qbd(J = 0), 
are obtained in combination with calculations in Ref. [lfj| (see 
text for discussion). All values are given in cm -1 . 



_J_ 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 



X(J) 



X, 



rol+QED 



(J) 



118.48684(10)^ 
354.3733(2)^, ~ 
705.5189(3£ 
1 168.7982(2^, 
1740.1895(3)^. 



0.00174(13) 

0.0049(2) 

0.0092(3) 

0.0157(2) 

0.0220(3) 



2 414 
3187 
4051 

5 001 

6 030 
7132 

8 298 

9 523 
10 801 
12123 
13 485 



898(5) 

472(5) 

943(5) 

963(5) 

921(5) 

066(5) 

600(5) 

.794(7) 

008(9) 

.83(1) 

99(1) 



0.031(5) 

0.040(5) 

0.049(5) 

0.058(5) 

0.069(5) 

0.081(5) 

0.087(5) 

0.101(7) 

0.103(9) 

0.12(1) 

0.13(1) 



a From Ref. [27J] , obtained from a fit of IR quadrupole transition 

frequencies. 
b Derived from EF — X transition frequencies from Refs. l2fj|422fl . 



ergy, signifying the ortho-para splitting quoted from Jen- 
nings et al. [27], to complete the rotational sequence. 

Relativistic and radiative corrections ^rei+QED(<^) are 
extracted from the difference in the experimental rota- 
tional level energies X(J) and those from the nonrela- 
tivistic calculations of Pachucki and Komasa [16]. We 
note that the A'icI+qed^) values, also listed in Table ITT1 
are taken with respect to X 1c i + qeb(J — 0). The value 
-X'roi+QED {J = 0) for the J = quantum state is equiv- 
alent (in magnitude) to the correction for the H2 dis- 
sociation energy calculated by Piszczatowski et al. [l4| 
to be +0.7283(10) cm^ 1 . Hence, the present investi- 
gation probes differential relativistic and radiative ef- 
fects in the rotational energy sequence. The value for 
X ic \+qed{J = 0) from Ref. 14] comprises a full calcula- 
tion of the leading-order relativistic correction and radia- 
tive corrections up to 0(a 3 ) orders in atomic units, and 
further includes estimates of radiative corrections of or- 
der 0{a A ) and recoil corrections up to 0(—a 3 ). These 

result in a value of 0.5319(3) cm -1 for the leading-order 
relativistic correction and 0.1964(9) cm -1 for radiative 
corrections (including higher-order relativistic terms) for 
X ic \+qed(J = 0). In the following discussions, we refer 
to the combination of relativistic and radiative correc- 
tions as QED corrections. The progression of rotational 
level energies X( J) and the experimentally-derived QED 
corrections AT rc i+QED(^) are plotted in Fig. O The QED 
corrections are in the order of 10 -5 smaller compared 
to the level energies, therefore, high experimental accu- 
racy is necessary to observe the QED effects. In addi- 
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FIG. 2. (Color online) Rotational level energies X(J) (tri- 
angles) of the H2 vibronic ground state X, v = 0. The rel- 
ativistic and QED corrections X tc \ + qw(J) (circles) are de- 
rived from the experimental X(J) and the calculated values 
in Ref. 16]. The unfilled squares (also in the inset) are from 
the calculations of X re l+QED{J) by Komasa et al. [17l | . 



tion, since X tc i + qed(J) increases with increasing J (for 
J < 22), it is possible to observe greater (differential) rel- 
ativistic and radiative corrections as higher J quantum 
states are probed. 

In Fig. [2] the experimentally determined values for the 
QED corrections are compared with the recent and yet to 
be pub lished calculations of X re i+QED(»/) by Komasa et 
al. [l7( for quantum states up to J = 16. The agree- 
ment between experiment and theory for the QED cor- 
rections of X 1 Yj+v = 0, J ground state rotational lev- 
els in H2 is remarkably good: this is the key result of 
the present study. The level of agreement suggests that 
the systematic uncertainty contribution to the transition 
frequencies, due to chirp phenomena in the PDA laser 
system, is overestimated in the present wavelength range 
used. The calculated QED corrections [with respect to 
-^rol+QED (J = 0)] for the full set of rotational levels (up 
to J = 31) in the X, v — band is plotted in the inset 

of Fig. m 

The above conclusion on the QED tests in H2 is in prin- 
ciple dependent on the correctness of the determination 
of EF x £+ ,v = 0, J level energies in the work by Bailly et 
al. [22| . A further test on the QED contributions relies on 
combination differences in the ground state, which we de- 
fine as the level splittings between J and J + 2 quantum 
states: A(J -> J + 2) = X(J + 2) - X(J). By com- 
paring EF — X (0, 0) transition frequencies belonging to 
the O(J), Q(J), and S(J) branches, as listed in Table Q] 
values for A( J —> J + 2) can be deduced in a straightfor- 
ward manner. Again, in order to address the pure QED 
effects A re i + QED(</ — > J + 2) in the molecule, the nonrel- 
ativistic contributions [16J to the combination differences 
A(J — ¥ J + 2) are subtracted. The resulting QED contri- 
butions A rc i + QED(</ — > J+2) are finally plotted in Fig. EI 
The upright triangles are derived from combination dif- 
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FIG. 3. (Color online) QED test for the combination dif- 
ferences in H2 level energies. Upright triangles: QED con- 
tribution to combination between O and Q lines; inverted 
triangles: QED contribution to combination between S and 
Q lines. The horizontal axis of the upright and inverted tri- 
angles are shifted slightly for clarity. The unfilled squares are 
from the calculations of Komasa et al. [17I ] . 



ferences between O and Q lines; the inverted triangles are 
derived from differences between S and Q lines. These 
quantities do not in any way depend on actual values of 
EF(J) excited state energies. It is noted that transi- 
tions connect only even-J levels or only odd- J levels in 
H2, as is always the case in para- and orthohydrogen. 
The corresponding theoretical values for the relativistic 
and radiative corrections A re i + QED(>/ — > </ + 2) for the 
same combination differences are also plotted in Fig. [3] 
for comparison. These purely theoretical QED contribu- 
tions were derived from combining the calculations for 
level energies in Refs. lg, |l7|. Figure [3] can be inter- 



preted as the first derivative, with respect to J, of the 
^ic1+qed(>/) sequence in Fig. ® thus explaining why the 
maxima in Figs. [2] and [3] are located at different J values. 
This second test of QED, based on combination dif- 
ferences between J and J + 2 levels, again shows perfect 
agreement at the level of estimated uncertainties. While 
the advantage of this second QED test is that it does 
not rely on previously measured EF(J) level energies, it 
has the disadvantages that significant cancellation of rela- 
tivistic and radiative contributions lead to smaller energy 
differences, and that the uncertainties in the experimen- 
tal data are slightly larger since combined errors must be 
taken. 

With the present study, the field of molecular spec- 
troscopy is opened up to include effects of quantum 
electrodynamics, i.e., the calculation of self-energy, vac- 
uum polarization, and high-order relativistic and radia- 
tive corrections is required for an accurate representation 
of level energies in molecules. A rotationally-hot popu- 
lation distribution of H2 molecules is created for a sys- 



tematic and precise spectroscopic study on the sequence 
of rotational level energies in the H2 ground electronic 
state. At the accuracy limit of the experiment (0.005 
cm -1 ) perfect agreement is found between experiment 
and theory on the QED contributions to level energies 
up to quantum state J = 16. Since the accuracy level of 
QED calculations is claimed to be an order of magnitude 
more precise than the present experiments [171 ] , there re- 
mains room for improvement to conduct more stringent 
experimental QED tests in H2. 
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